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O ; Abstract 

The dynamics of physical theories is usually described by differential equations. Difference 
equations then appear mainly as an approximation which can be used for a numerical analysis. As 
such, they have to fulfill certain conditions to ensure that the numerical solutions can reliably be 
used as approximations to solutions of the differential equation. There are, however, also systems 
where a difference equation is deemed to be fundamental, mainly in the context of quantum gravity. 
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. Since difference equations in general are harder to solve analytically than differential equations, it 
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can be helpful to introduce an approximating differential equation as a continuum approximation. 
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loop quantum cosmology, we show that the conditions are restrictive and serve as a selection 
criterion among possible quantization choices. 

PACS numbers: 0460P, 0460K, 9880Q 



In this paper implications of this change in view point are analyzed to derive the conditions 
that the difference equation should satisfy. The difference equation in such a situation cannot be 
chosen freely but must be derived from a fundamental theory. Thus, the conditions for a discrete 
formulation can be translated into conditions for acceptable quantizations. In the main example, 
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I. INTRODUCTION 



While differential equations play an outstanding role throughout physics, it has long 
been expected that a quantization of gravitational systems will lead to a discrete structure 
of space-time. Then, the differential equations must be replaced by difference equations at 
a fundamental level. The process of discretization itself is, in fact, well-known since it is 
used to study differential equations numerically. In this case the difference equation is only 
an approximation to the differential equation which is better suited to a numerical analysis. 
When the fundamental equation is discrete, however, the point of view must change. Now, 
analytic calculations are more complicated and less exact techniques are available than would 
be for a differential equation. Thus, it becomes convenient to consider, in certain regimes, a 
differential equation as an approximation to the difference equation. This is also necessary 
since the continuum formulation gives an excellent description of the real world such that 
there must be a continuum approximation to the discrete theory. 

In both cases the relation between a difference and a differential equation is analyzed. 
But there is an important difference: In numerics Q one has a large freedom in choosing 
the difference equation which can be exploited to obtain a very good approximation to the 
differential equation with helpful properties for the numerical analysis. If the difference 
equation is fundamental, e.g. if it follows from a basic quantum theory, then it cannot 
be changed so freely. Only quantization choices and ambiguities can be exploited which 
are often more limited and have to respect other consistency conditions of the theory. In 
particular, while the order of the difference equation should not be less than that of the 
differential equation to allow the same freedom as in the continuum formulation, there is 
no guarantee that it is not of considerably higher order. Also in numerics higher order 
schemes are being used to improve the numerical results, but one has to choose the higher 
order discretization very carefully. Whenever one uses a higher order, there are additional 
solutions, unwanted from the numerical point of view, which have to remain suppressed 
throughout the evolution. When the fundamental difference equation is of higher order than 
a continuum differential equation, there will also be additional solutions which correspond 
to a new freedom of the theory not seen in the continuum picture. While those solutions can 
play a physical role in certain regimes, their role must be limited in semiclassical regimes 
for the continuum formulation to be valid. Now, one has to check if the freedom allowed 
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by different quantization choices can not only provide a discrete equation which has the 
expected continuum approximation, but also, in the case of higher order, does not lead to 
growing modes which would soon dominate the continuum approximations. 

The aim of this paper is to extract conditions for a fundamental difference equation to 
have a good continuum picture in this sense. Moreover, we will discuss how these conditions 
can be used as selection criteria among different quantizations or ambiguities. It turns out 
that they can be quite restrictive, thus having implications for quantization issues which at 
first glance seem only indirectly related to the fundamental law. 

Our motivation, and the main example for a theory where in fact a discrete fundamental 
evolution law has been derived, is loop quantum cosmology. The issue of discretization 
itself as a consequence of quantum theory or gravity has been studied before [2|, but in 
particular the issue of higher order has not always been appreciated. In order to cover 
also other examples of a discrete formulation, we will hold the discussion more general, but 
will use the example of loop quantum cosmology as an illustration throughout the paper. 
Moreover, we expect that loop quantum cosmology, or more generally quantum geometry 
which it is a part of, will be the prime application of the selection criteria. The reason 
is that quantum geometry is the most advanced formulation which has led to a discrete 
structure. Its kinematical setup is at a mathematically precise and well-defined level which 
demonstrates the main quantization ambiguities involved in the definition of the basic laws, 
mainly the Hamiltonian constraint. In such a scenario the selection criteria will be most 
useful since they can be tested explicitly and can restrict the possible choices. While explicit 
calculations in the full theory are still outside current reach, one can study the basic law 
explicitly in homogeneous models. Most of them, those with non-zero intrinsic curvature, 
require an additional input compared to the full theory such that the contact is not so close 
and more ambiguities arise [3j. It is then helpful to use the criteria discussed here in order 
to select a quantization. 

In section II we specify the class of difference equations considered. We briefly recall 
the quantization of homogeneous models and describe the form of the ensuing fundamental 
difference equations. 

In section III, we discuss the intuitive notion of continuum approximation and its for- 
mulation in terms of pre-classical solutions. For definiteness and clarity, in this section we 
discuss fundamental equations which are ordinary difference equations. The notion of pre- 
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classicality is carefully defined without appealing to any limiting procedure (Def. |2J). Its 
implications for the difference equation are detailed. Subsequently a constructive procedure 
is given to demonstrate the existence of pre-classical solutions. This procedure naturally 
leads to a notion of local stability which is defined here (Def. |5J). The conditions of ad- 
missibility of a pre-classical solution together with local stability are given in terms of the 
coefficients of the fundamental difference equation. These are our main results summarized 
in Props. ^ and El 

In section IV, the case of a fundamental partial difference equation is analyzed. The 
analysis closely follows the steps detailed in section III. In some cases the fundamental 
equation can be interpreted as a non-homogeneous ordinary difference equation and this 
can be done in several ways. For each of these, the definition of local stability from section 
III can be applied. The fundamental difference equation can then be constrained by requiring 
that the local stability criteria must be satisfied for any evolutionary interpretation. The 
example of homogeneous loop quantum cosmology p is described briefly. 

In the concluding section V, we summarize our results. 

An appendix is included which illustrates an example of an older, unstable quantization 
of isotropic loop quantum cosmology. 

II. DIFFERENCE EQUATIONS 

In this section we fix our general notation of the class of difference equations considered, 
and describe the more special way in which they arise from loop quantum cosmology. 

A. General set-up 

We assume that there is a continuum formulation for the system under consideration 
which has a number of continuous parameters p 1 (which can be space or space-time coordi- 
nates, or in the case of cosmology minisuperspace coordinates) and a linear, homogeneous 
differential equation T>s = of order k for a (wave) function s(p 7 ). These basic variables 
have an underlying discrete structure given by the lattice 




(1) 
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where m; are integers and 7 is the dimensionless parameter which controls the scale of the 
discreteness. The parameter t will always be fixed and can absorb possible dimensions. At 
the discrete level, the fundamental equation is a difference equation of the form 

k 

E <; } (7) = (2) 
{ii=~k} 

of order 2k for a (wave) function S{ mj } living on the lattice (the notation {. . .} in the 
sum indicates that all ij are summed over). The coefficients can in general depend on the 
coordinates mj as well as the parameter 7. 

Due to the relation between p 1 and m I we expect a relation between the discrete and the 
continuous formulation when 7 is small and the mi are large, which will be assumed from 
now on. If 7 is small but non-zero, the differential equation can only be an approximation to 
the difference equation; a strict relation could be achieved in the limit 7^0, mj — > 00 such 
that p 1 is unchanged, which we will call, as in loop quantum cosmology, the pre-classical 
limit. Since a relation between the formulations will be demonstrated by using expansions 

S{ mi +ij} = s{p I {m I )) + S ^2iiid I s(p I (m I )) H (3) 

one also has to assume properties of the wave function concerning its oscillation on small 
scales. This will be discussed in more detail later. 

We will also assume that the coefficients of the difference equation can be expanded in 
terms of 7, 

00 

<} } (7) = E^4!}(p'N) (4) 
j=-i 

where poles at most of order I appear. (Some, but not all of the l \ (p 1 ) can be zero.) Note 
that we use the p 1 to parameterize the coefficients in the expansion. This is necessary because 
later we will use the expansion in the continuum limit where some of the 7-dependence is 
absorbed in the toj. Furthermore, we assume the following conditions for the coefficients to 
be satisfied at least for large mi which would be used in a pre-classical limit or approximation 
(at small raj, the conditions can be violated by quantum effects): 

unitarity conditions A^f = (^{^ } } )* (5) 
parity conditions = ^^{i 7^ ■ (^) 
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The first condition ensures that the equation with a reversal of the direction of "time" is 
solved by the complex conjugate of the original wave function: 

k k / k \ * 

A {U) S \"H-il} = Z2 A {-ii} S {mi+ii} = I A {h} S (™/+ i /} = 

{il=-k} {h = -k} \{i I = -k} J 

while the second one guarantees that a solution remains a solution when all the coordinates 
are reversed. 

Depending on the situation and the physical meaning of the parameters p 1 , we will usually 
assume rotational symmetry of the equations, i.e. permutation symmetry for the coefficients 

771/ . 



B. Loop quantum cosmology 
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In diagonal homogeneous loop quantum cosmology we have three independent vari- 
ables p 1 , I = 1, ... ,3, which are the components of the densitized triad. They determine 
the spatial metric as ds 2 = a 2 - (a/) 2 where u 1 are fixed invariant 1-forms which specify the 



particular model, and aj = y\p J p K ' jp I \ if £ijk = 1- If we restrict to isotropic models, there 
is a single p and a = \f\p\ i s the scale factor. The discreteness is controlled by the Barbero- 
Immirzi parameter 7 of quantum geometry, and £ = \t\ with the Planck length £p (or 
£ = \i\ in the isotropic case). Note that the regular lattice discretization p I {mi) = ipfi%mi 
is not an assumption but derived from the theory via the spectra of densitized triad oper- 
ators; it is a general feature of quantum geometry. For diagonal homogeneous models the 
lattice structure follows from the eigenvalues p I {m{) of triad operators p 1 with eigenstates 
|mi, 777,2,7713). A wave function is supported on this lattice, s miim2im3 in the triad represen- 
tation defined by \s) = J2{ mi } s m 1 ,m 2 ,m 3 \ m i, m 2, m 3) for an arbitrary state \s). Since the 
states |mi,m 2 ,m 3 ), defined as eigenstates of the triad operators, have the freedom of an 
arbitrary (possibly m/-dependent) phase factor, the triad representation and the fundamen- 
tal equation for a wave function are fixed only after a choice of the phase factors has been 
made. 

Another essential part of quantum geometry is that the densitized triad is conjugate 
to a connection, which in the diagonal homogeneous case also has three components cj = 
Tj — jKj where Tj is the spin connection and Kj the extrinsic curvature. It is quantized in 
quantum geometry not directly, but via its holonomies P|. In the homogeneous case, these 



holonomies are given essentially by sin(|cj) and cos(|c/) acting as multiplication operators 

in the connection representation. If we choose the phases for the triad eigenvectors as 

3exp(i(miCi + m 2 c 2 + m 3 c 3 )/2) 
(c 1) c 2 ,C3|m 1 ,m 2 ,m 3 ) = 2 2 sin ( Cl/2 ) sin ( C2 / 2 ) sin ( C3 / 2 ) ^ 

the holonomy operators act on a wave function S{ m/ } via, e.g., 

(sin( 2^l)^)mi,Ti2,r7i3 2 ^("Smi+l,r7i2,rri3 S mi _ 1,7712,7713) (8) 

(cOs( 2Cl)'S) mijm2j rri3 2 (^ m i+li m 2i m 3 ^7711—1,7712,7713) • (9) 

The fundamental equation is given by the Hamiltonian constraint equation for a wave 
function. In the full theory, a part of it is constructed with a Wilson loop evaluated in 
the connection. It can straightforwardly be specialized to a homogeneous model provided 
that there is no intrinsic curvature (i.e., the Bianchi I model). To respect the symmetry 
we are only allowed to use integral curves of vector fields generating the symmetry (dual 
to the invariant 1-forms used above) in forming the Wilson loop; this can be done by 
forming a square from four such curves which close because the vector fields commute. 
The Hamiltonian constraint equation for the Bianchi I model is then essentially fixed taking 
the form (fTT?|) with non-zero coefficients 

A™ e \'™el'fl 3 (l) — e l e 2 7 3fi 1 £p 2 (Vm 1 -2e 1 ,m 2 -2e2,m :i +l — Vm 1 -2e 1 ,m 2 -2£ 2 ,m a -l) (10) 

for ej = ±1, and similarly for other coefficients by permutation symmetry. We used the 
volume eigenvalues 

V muma ,ma = (lT4)^Vl^l^2W3| (11) 

to simplify the notation, £p = kH, and k = 8ttG is the gravitational constant. If there is 
matter or a cosmo logical constant, there will also be a coefficient (it can be operator valued, 
acting on the matter dependence of the wave function) 

A™fr™ 2 ' m3 (<y) = //matter (mi, m 2 , m 3 ) • (12) 

In the notation of (|1J we have a pole of order I = 2 and non-vanishing coefficients which can 
be read off from 

Kin- 2ei,m2— 2e2,T7i3+l •mi- 2ei, 7712— 2e2, 7713— 1 

| 7 4V|pV/p 3 | (1 - WUei/p 1 + e 2 /p 2 ) + 0( 7 2 )) . 
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The order of the partial difference equation is 2k = 4. 

The fact that we have a difference equation can be traced back to the Wilson loop 
operator which can be decomposed into a product of the operators JSJ), Q giving us the 
basic difference and mean. It also demonstrates the point that a fundamental difference 
equation can easily be of higher order: the basic difference ((SJ) is a second order difference 
operator, rather than first order (corresponding to a leap-frog scheme in numerics). It is a 
basic result of the theory and cannot be changed to simplify the equations. 

In other Bianchi models with non-zero intrinsic curvature the situation is not so clear-cut. 
Their symmetries are generated by vector fields which do not commute, thus not forming 
closed loops in the way described above. We have to correct for this fact by making another 
choice which cannot be motivated from the full theory. There needs to be an additional step 
because non-zero intrinsic curvature in a model is realized by a non-zero spin connection. 
While this is an invariant statement in the model, it does not make sense in the full theory 
where the spin connection can be made to vanish locally. Therefore, it does not pose a 
oroblem in the full theory, but we have to face it in any model where it does not vanish 
sj]. Guidance by the full theory, which we could use for the Bianchi I model, is lost which 
results in more quantization choices and ambiguities. 

The purpose of the conditions studied in this paper is two-fold. We have to check if they 
are fulfilled in the Bianchi I model which then, due to its close relation to the full theory, can 
be regarded as a consistency check for the full quantization. For the other models we can 
use the conditions to select among the possibilities for a quantization (if they are possible 
to realize at all). 



C. Pre-classical approximation vs. pre-classical limit 

It is obvious from the basic definitions that the smaller the parameter 7, the closer we 
are to a continuum formulation. This also requires to have large values of the discrete 
variables mj, and only small variations of the wave function S{ m/ } between adjacent lattice 
points. If we perform the continuum limit 7^0 exactly, then we must simultaneously 
lave mi — > 00 and s must become a smooth function s(p 7 ); this is the pre-classical limit 
(|. From a different conceptual point of view, we can also work with small but non-zero 
values of 7, which is the physically correct way in the case of loop quantum cosmology 
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where 7 = loe,(2)/7r\/3 can be inferred from calculations of black hole entropy in quantum 
geometry In this case we still have to consider large mj and small oscillations of the 

wave functions at Planck scales for the continuum •pre- classical approximation to be valid. 
This corresponds to a regime of large volume and small curvature in the case of gravity. 

Concerning the main theme of this paper, the first point of view of a pre-classical limit is 
the one usually taken in numerics. When the differential equation is regarded as fundamental 
and the difference equation as an approximation, the limit has to be performed eventually 
such that 7 can be assumed to be as small as necessary for a particular mathematical 
result. Thus, mathematical techniques developed for the numerical analysis of differential 
equations can be used when studying the pre-classical limit. When the difference equation is 
fundamental, however, 7 has a physical value which is non-zero. In such a situation results 
which rely on the fact that 7 can be chosen as small as appropriate, will not apply and new 
definitions and theorems have to be developed. In this paper we present a first approach 
giving basic definitions and some implications. 

III. CONDITIONS ON THE FUNDAMENTAL DIFFERENCE EQUATION 

To simplify the presentation let us first work in the context of isotropic loop quantum 
cosmology l| so that we have a single triad component p with eigenvalues p(m) = ^ipm 
and eigenstates \m). The fundamental difference equation is of the general form (J2J) with 
only one summation label, 

A™( 7 ) Sm+fc + A^_ x {i)8 m+k . x + ■■■ + A™ k ( 7 )s m _ fe = , (13) 

where the coefficients satisfy the unitarity conditions for large m, 

A-i * (AT)* . (14) 

The 7-expansion of the coefficients is such that only one pole term of order I = 2 appears, 
while all other terms have non-negative power of 7. 

Note: In practice, the difference equation resulting from a quantization of the Hamiltonian 
constraint is not unique. As discussed in the preceding section, the discrete wave functions s m 
are defined only modulo m-dependent phase transformations which also alters the coefficients 
of the difference equation. Such a freedom has to be kept in mind while checking the 
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conditions we will derive for any particular difference equation. For homogeneous difference 
equations which arise in the isotropic context, the coefficients are defined only up to common 
factors. Additionally, for large m in particular, such common factors do occur but can be 
canceled. In the following, the coefficients will be assumed to have no such common factors. 

A. Continuum Approximation 

When a discrete formulation is used at the fundamental level, it has to reproduce the 
familiar continuum idealization at large scales or in semiclassical regimes as a very good ap- 
proximation. This presents a necessary condition for any discrete formulation; in particular 
a quantization of geometry has to reproduce the familiar continuum metric formulation as 
an approximation when quantum effects can be expected to be negligible. The latter is real- 
ized when the length scale on which the continuum formulation is used is large compared to 
the Planck scale on which quantum geometry is to be used, and curvature is not significant. 
In the context of homogeneous cosmology the corresponding region of phase space has large 
values of the triads (and by implication of volume) and small values of the extrinsic curva- 
ture. The intrinsic curvature of a minisuperspace model, by contrast, is a fixed function on 
the phase space and thus cannot be independently required to be small in a semiclassical 
regime; it can be large even in regions where the volume is large and the extrinsic curvature 
is small, which has to be taken care of properly. 

For gravity, the familiar Wheeler-DeWitt quantization of homogeneous models jloL lri| 
is obtained using a Schrodinger quantization of the continuous geometry variables, i.e. the 
spatial metric and the extrinsic curvature. In this quantization, the fundamental formulation 
is still continuous: wave functions are functions of continuous labels p and the constraint 
equation leads to the Wheeler-DeWitt differential equation. Since the Wheeler-DeWitt 
quantization already uses continuous geometric variables, it is convenient (but not necessary 
since, after all, the Wheeler-DeWitt quantization is not well-defined for the full theory) to 
require that the wave functions of this quantization together with the Wheeler-DeWitt 
equation they satisfy should be reproduced as a suitable limit or at least as the leading 
behavior in an approximation from some solution(s) of the fundamental difference equation. 

The indication of the regime of semiclassical behavior and the justification for matching 
with the Wheeler-DeWitt quantization can be further seen by observing that in a triad 
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representation, the basic difference operators of the discrete formulation are given by sin(|c) 
as in (|SJ), while the basic differential operators of a Wheeler-DeWitt quantization are given 
by c = i'jipd/dp. Thus, if c can be assumed to be small, we in fact have sin(^c) ~ \c and the 
difference operators immediately yield differential operators for small c. This assumption, 
however, is justified only in the simplest cosmological models, the Bianchi I model and its 
isotropic sub-models in which the intrinsic curvature, represented by the spin connection T, 
is zero and c is essentially the extrinsic curvature. Otherwise, c also has a contribution from 
intrinsic curvature and the relation to the continuum formulation is less direct; its proper 



treatment is discussed 



in [i, 



121. 



The condition of large volume is easier to implement since the large volume regime nat- 
urally corresponds to large eigenvalues of the triad operator, m 3> 1. If finite increments 
of the label m appearing in the difference equation have only small effects on the wave 
function (which also depends on the curvature) then one can expect difference operators 
acting on sufficiently slowly varying functions of m to be well approximated by differential 
operators. In this way, one can conceive the possibility of recovering the Wheeler-DeWitt 
quantization. Thus, the basic idea is to see if there are 'sufficiently slowly varying' discrete 
wave functions, s m which can be interpolated by 'sufficiently slowly varying' solutions of 
the Wheeler-DeWitt equation. The requirement that the Wheeler-DeWitt quantization be 
obtained as a good approximation translates into the requirement that the fundamental 
difference equation admit at least one such slowly varying solution. We will now make these 
statements precise and deduce their implications. 



1. Preliminaries 



In the following, the domain of validity of a continuum approximation is specified by 
1 <C M m i n <C M max and m e (M m i n ,M max ). For mo in the domain of validity, we will 
specify its Am- neighborhood as the interval [mo, mo + Am]. Statements valid over Am- 
neighborhoods will be termed as local statements. Typically, 1 <C Am <C m will be 
assumed. In some situations the domain of validity can extend over an infinite range of m. 

The coefficients in the difference equation depend on 7 and m. The 7 dependence is such 
that A™ = 7- 2 (A(" 2 ))™ + (A^)™ + 7(AW)™ + ■■■. Note that there is a term of order 7" 2 
while there is no 7 _1 term in the coefficients, as discussed in Section Til Bl All the (A^)™ 
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also depend on m which, in a Am neighborhood of some mo in the domain of validity, can 



be expanded as: (A^)™ 0+Sm ~ + ^Bf ] + ■■■ , < 5m < Am , where Af\ B\ r ] are 



constants depending on the fixed mo but independent of 5m. 

A few preliminary definitions: A function ip of a continuous variable p is said to be locally 
slowly varying around po if ip{p) can be Taylor expanded about po with succeeding terms 
smaller than the preceding terms for p G [po,Po + Ap]- A bound on the ratios can be used 
to quantify degree of slow variation but we will not need it. We will typically need terms 
only up to second order. Correspondingly, a sequence s m is said to be locally slowly varying 
around mo if it can be obtained via s m := ip{j>{m)) where p{m) = ^ji^m, ip(p) is locally 
slowly varying. In other words, a sequence is slowly varying if it can be interpolated by a 
slowly varying function, which, of course, will not be unique. 

Notice that the definition of the variable p(m) contains an explicit 7 dependence. Thus, 
even if a slowly varying ip(p) has no 7 dependence, the corresponding slowly varying se- 
quence does depend on 7 and is in fact Taylor expandable in 7. Conversely, if a slowly 
varying sequence has a specific type of 7 dependence (because it is a solution of the differ- 
ence equation containing 7, for instance), the corresponding interpolating function can be 
independent of 7. 

Thus, a slowly varying function and a slowly varying sequence can be written as 



Let us assume now that the fundamental equation admits solutions which are locally 
slowly varying around some m in the domain of validity. 

2. Pre-classicality 

Using the 7-expansions of the coefficients and substituting (JTfij) for a locally slowly varying 
s m+i in the fundamental difference equation (|T3j) leads to, 

k 




(r 2 (A^)r + {a^)t + 704 (1) )r + • • •) 



i=—k 




(17) 
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It is clear that we can get a differential equation for ip{p) in this manner for every locally 
slowly varying solution s m of the fundamental equation. 

Essential for a continuum approximation is not just to get a differential equation but 
that a differential equation so obtained should be independent of 7. The 7 independence 
condition is necessary since by definition 7 controls the discreteness such that the continuum 
variables are obtained for 7 — > 0. For instance, in the case of gravity the differential equation 
is expected to be the standard Wheeler-DeWitt equation which knows nothing about 7. 

It follows immediately that ip{p) must satisfy equation(s) obtained by equating coefficients 
of powers of 7 to zero. There are, of course, infinitely many such equations. 

From any of these equations, one can derive a 'Hamilton- J acobi' equation for the phase 
of the wave function in the usual manner and read off a classical Hamiltonian constraint. 
Only one of these equations reproduces the correct classical Hamiltonian constraint. In view 
of the small value of 7, this equation satisfied by ip{p) is to be considered as the leading 
approximation and consequently, the equations corresponding to lower powers of 7 must be 
identically zero without any conditions on ip{p). 

Thus we define: 

Definition 1 A locally slowly varying solution around m of the fundamental difference 
equation is said to be locally pre-classical around m Q if there is an interpolating slowly 
varying function ip{p) which solves a differential equation that has no dependence on 7 and 
has the correct classical limit. 

Thus, the 7 independence condition together with the correct classical limit, implies that 
the o(7~ 2 ), o(7 _1 ) terms must be zero identically and the 0(7°) term must contain derivatives 
of ip(jp). Furthermore, in a Am-neighborhood of mo we can use (A^)™ ~ A ( [ ] which also 
removes the explicit m dependence. This leads to: 

k k k 

A { r 2) = 0=J2 iA< f 2) > and E A i~ 2) i2 ^° v m e ( m °' m ° + Am ) ' ( 18 ) 

i=—k i=—k i=—k 

and we obtain, to the leading order in 7 and the approximating differential equation as 
14 A ( " 2) ^ + A (0) ) m + o(7, Am/m ) = . (19) 
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The approximating differential equation will precisely match with the Wheeler-DeWitt 
equation if the ratio of the coefficients of the two terms on the left hand side of (fTU|) matches 
with those of the Wheeler-DeWitt equation. This indeed turns out to be the case. We 
will refer to the leading order approximating differential equation as the Wheeler-DeWitt 
equation. Using the Wheeler-DeWitt equation for capturing the continuum behavior is 
thus an approximation. From the perspective of the discrete solutions, locally pre-classical 
solutions are still exact solutions of the fundamental equation. 

Thus we have seen that the admissibility of even a single locally pre-classical solution 
puts strong restrictions on the coefficients of the fundamental equation: 

Proposition 1 If a difference equation of the form (Q3J) with coefficients having a second 
order pole in the parameter 7 admits a locally pre-classical solution, then the leading order 
terms of the coefficients, Al~ 2 \ must satisfy (03). 

We are now ready to define pre-classical solutions. 

Definition 2 (Pre-Classicality) A solution of the fundamental equation is said to be pre- 
classical if it is locally pre-classical around every m in the domain of validity. 

The requirement of the existence of a pre-classical solution implies that all the local 
statements above, in particular the conditions in (J18|) must hold throughout the domain of 
validity. 

Notice that these conditions are homogeneous in the coefficients. For the equations of the 

( — 2) 

isotropic homogeneous models, there are only three non-zero coefficients A\ and the two 
conditions of (|18|) fix these up to an over-all scaling. It is the over-all scaling that contains 
the m dependence. The remaining ratios being independent of mo are the same for local 
neighborhoods around any mo in the domain of validity. 

To summarize: We have reformulated the definition of a pre-classical solution without ref- 
erence to any p r .cla,ical Bndt Q. It is m0r e realistic s,„ C e it worfcs with Sm a„ bu t _ 
7 and large but not infinite m, captures the essential ideas of getting a continuum descrip- 
tion as an approximation in terms of solutions of the Wheeler-DeWitt equation and directly 
gives restrictive conditions (|18|) on the coefficients of the fundamental difference equation. 
One can construct approximate pre-classical solutions from (slowly varying) solutions of the 
Wheeler-DeWitt equation as outlined in 

14 



B. Local Stability 



We now address the issues of the existence of pre-classical solutions, the number of inde- 
pendent pre-classical solutions and the construction of approximate pre-classical solutions. 
We will work exclusively within the discrete formulation and assume that the coefficients of 
the fundamental equation satisfy the conditions of (|18|) which are necessary conditions for 
the existence of a pre-classical solution. As a consequence, we will construct local approxi- 
mations and use them to argue about the existence of pre-classical solutions. 

1. Approximate pre-classicality 

Consider locally slowly varying solutions (pre-classical ones being a special case) of the 
fundamental equation. We know that such solutions are Taylor expandable in 7 : s m = 
Sm + 7 Sm + 7 2 Sm + • • •■ Substituting in (JT3J) and using the fact that the coefficients are 
locally almost constant gives, 

= 7" 2 (t + T' (E ^<&) 

\i=-k J \i=-k J 

\i=—k i=—k / 

It is evident now that a locally pre-classical solution can be constructed as a power series in 
7 by choosing Sm , Sm to be solutions of 

k 

Y,A ( f 2) s m+t = 0, (21) 

i=—k 

(2) 

while Sm is a solution of the non-homogeneous equation 

E 4- 2 '^, + E A f S m+i = • (22) 
i=—k i=—k 

The higher order terms are to be constructed recursively from solutions of non-homogeneous 
difference equations with constant coefficients. Truncating the Taylor series to any finite 
order constitutes an approximation. Thus we define: 

Definition 3 A solution of the form s m = Sm + 7 Sm + 7 2 Sm with Sm satisfying \21^2ty 
is said to be a 2nd order approximate, locally pre-classical solution or 2nd order local ap- 
proximation for short. 
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Apart from being Taylor expandable in 7, locally slowly varying sequences also have the 
form s m+ i ~ a m + ib m + i 2 c m + ■ ■ •, where b m ~ 0(7), c m ~ o(7 2 ) etc. as seen in (|TB|) . In fact, 
the equations (|2T1 do admit solutions of this form. To see this let us take the ansatz 
(for r = 0, 1, 2 ) , 

s^ o+6m ~ a (r) (m ) + (m + 5m)b {r) (m ) + (m + 5m) 2 c {r) (m ) H , < 5m < Am . (23) 

Substitution in the equations (|2*T| E2J) implies that a^ r \m ), b^(m ) are arbitrary, 
c' r '(mo) = for r = 0, 1 and c' 2 '(mo) is determined in terms of the a'°'(mo ), &(°)(m ) 
and coefficients of the difference equations. The conditions ()18|) are crucial for these 
statements. We can choose the a^(mo),l)' r '(mo) coefficients so that s mo+ $ m is obtained 
as a Taylor series in 'jSm. Explicitly the choices are: b^ (mo) = 0, a^(mo) = 
—rriob^ (mo) , a' 2 '(mo) = — mflc' 2 '(m ), fe( 2 )(m ) = —2m c^(m ). This leaves us with 
two free parameters, a'°'(mo), M^mo). We have thus constructed 2nd order local approx- 
imations (in a Am neighborhood of mo and to order 7 2 ) which are parameterized by two 
constants and given by, 

s mo +5m = a (0) (m ) + 7^& (1) (mo) + l 2 Sm 2 c {2) (m ) . (24) 

The pre-classicality properties of these solutions are obvious. It is also clear that the 
number of independent, exact pre-classical solutions can not be more than two since their 
local approximations are parameterized by two parameters (this is true even if we go to 
higher orders in 7). 

We can construct such two parameter families of approximate, locally pre-classical solu- 
tion around various rrij in the domain of validity. These should somehow be tied together to 
constitute an approximation to a single (common) exact pre-classical solution, as discussed 
in what follows. 

2. Construction of approximate, pre-classical solutions 

Now we give a procedure to construct a solution which will approximate a pre-classical 
solution in various local neighborhoods. To do so, we divide the domain of validity into 
several, non-overlapping local neighborhoods around m := M min ,mi, ■ ■ ■ ,m N := M max 

(—2) 

within each of which we have the constant coefficients A] (mj) satisfying the conditions 
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(|18|). We can construct local approximations in each of these local neighborhoods as outlined 
above. These local approximations now have to be 'matched' across the boundaries of the 
adjacent neighborhoods. 

Both the exact (|13j) and the approximated (J211 evolution equations are of order 2k, 
and to determine any particular solution we need 2k consecutive s m 's to be specified as 
initial conditions. We will generically refer to the set {s m , s m+ i, • • • , s m+ 2k-i} as initial data 
at m. In each neighborhood, the local approximations induce initial data of the form 

s mj+i = a (0) K) + lib^imj) + 7 Vc (2) (m j ) + • • • , i = 0, 1, • • • , 2k - 1 . (25) 

Likewise, we also have the 2k consecutive values from s^.+Am— 2k as > 

s mj +A mj -2k+i = a {0 \m j )+^(Am j -2k + i)b {1) (m j )+-f 2 (Am j -2k + i) 2 c {2 \m j ) + --- . (26) 

The initial data at rrij + i = rrij + Amj however can not be arbitrary since the exact 
equation holds across the neighborhoods. We can still use the constant coefficients to connect 
these values, however the constant coefficients will now involve a mixture of coefficients from 
both the adjacent neighborhoods and this mixture need not satisfy the conditions (|18|). 
Equivalently, the initial data so constructed at rrij + Arrij need not equal the data induced 
from the local approximation in the neighborhood [m J+ i,mj + i + Auij + i\. We will refer to 
the former as the extended initial data, the latter as the local data and the corresponding 
solutions of (|2"T1 l2~2"j) as extended solution and local solution. Note that the initial data consist 
of 2k values while we have a freedom of only two parameters in the local data. Since usually 
2k > 2, the two data sets are generically distinct and so are the two solutions. In particular, 
the extended solution cannot be a local continuum approximation. Using a local solution 
instead of an extended solution thus introduces errors which can accumulate with subsequent 
extensions to cover the full domain of validity. The error accumulation can even jeopardize 
the existence of a pre-classical solution itself. This is prevented by requiring 'stability' for 
the approximated evolution (J21] l22|) . 

To arrive at a suitable notion of stability for (|2"T]) . let us concentrate on the equation: 

k 

AiS m+i = , V m G [rrij + k, rrij + Am — k] (27) 

i=— k 

where A$ are constants which are assumed to satisfy the conditions (j!8)) . With the identifi- 
cation Ai := A\~ 2 \ this is the same as (}2"Tj) and has already been referred to as the locally 
approximated evolution. 
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It is convenient to denote the initial data by a vector in C 2fe , (S(mj)) 1 := s m . + j , i = 
0, l,---,2k — 1. On the space of initial data, we have the natural Euclidean norm 

c 2 M|5K-)|| 2 := E-=o 1 K-+il 2 - 

Consider now the approximated evolution of two 'nearby initial data'. Let S(rrij) be fixed 
initial data and consider other arbitrary initial data S'(mj) such that \\S'(mj) — S(mj)\ \ =: 
S mj < A. These data will generate two solutions s' m , s m whose difference, 6s m , will also 
satisfy the approximated evolution equation. We would like to see how much the two 
solutions differ at rrij + Am — 2k. 

The general solutions of difference equations of constant coefficients are given by linear 
combinations of elementary solutions of the form (z a ) m m Ta where z a is a root of the charac- 
teristic polynomial, ^ Aiz\ with multiplicity d a and < r a < d a — 1 Q]. Notice that (|T8|) 
implies that z = 1 is a root with multiplicity 2. Clearly, 5s m = ^2 a Y2 r a =0 SC ajra [z a ) ra m Ta 
where a ranges over the distinct characteristic roots. Since S'(rrij) is arbitrary, consider a 
5s m such that SC a ^ ra is non-zero only for a particular elementary solution. Then we have, 

S mj +Am-2k ^ \Za\ Am - 2k (l + r a Am ' 2k S mj . (28) 

Clearly, if there is a characteristic root with absolute value greater than one, then the 
two solutions deviate exponentially. Taking the data S m . to be local data and S' m . to be the 
extended data, we see that even if the two data are within a A neighborhood initially, they 
could evolve far apart if \z a \ > 1. This leads us to the definition: 

Definition 4 The approximated evolution \2T\ j is said to be stable if S mj+ i < A(l + 
o(^i)) V < % < Am - 2k for all 6 mj < A. 

A stable evolution then guarantees that the extended solution will be nearby the local 
solution i.e. \Ss m \ < A. It implies that all characteristic roots must have absolute value less 
than or equal to one. The coefficients of the difference equation also satisfy the complex 
conjugation property (JHJ) which implies that if z is a root, so is (z*)^ 1 . Combined with the 
stability requirement this implies 

Proposition 2 The approximated evolution of a difference equation (Q3J) with coefficients 
satisfying the unitarity conditions is stable if and only if all characteristic roots have 
absolute value one. 
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As noted above, in the 'first neighborhood' ([mo, mo + Amp]), we are left with a single 
free parameter, say a'°'(mo), due to the consistency condition. In the next neighborhood 
we have two free parameters, a'°'(mi), &W(mi). We can fix these by requiring that the 
extended and the local data at mi differ by the least amount eg. by a least square fit. This 
determines the free parameters and also determines S mi . We can take A = S mi . If the 
locally approximated evolution is stable in the second neighborhood, then the local solution 
in the second neighborhood is guaranteed to be within 8 mi of the solution extended from the 
first neighborhood. If locally approximated evolutions are stable for all the neighborhoods, 
then we can continue our procedure and construct a sequence s m with the property that 
it is an exact solution up to m and is an approximate, locally pre-classical solution in 
every neighborhood of the chosen partition of the domain of validity. It is thus essential 
for the construction of an approximate, piece-wise, pre-classical solution that all the locally 
approximated evolutions be stable. This leads us to define: 

Definition 5 (Local Stability) The fundamental difference equation is said to be locally 
stable if the locally approximated evolution in any local neighborhood around any tuq in the 
domain of validity, is stable. In particular this implies that the coefficients of the funda- 
mental equation satisfying ilfy) have to satisfy the further condition that their corresponding 
characteristic polynomials have all roots with absolute value one. 

In other words, 

k 

AS" 2) (m )^(m ) = implies \z(m )\ = 1, V m G [M min ,M max ]. (29) 

i=— k 

Note that the construction of an approximate pre-classical solution depends on the choice 
of M m in (up to which exact evolution is used) and also on the partitioning of the domain 
of validity. Hence such approximate pre-classical solutions are not unique. This is also true 
of exact pre-classical solutions. Since one does not know what initial data at m = would 
characterize a pre-classical solution, we can choose the local data at any mo G [M m j n , M max ] 
and evolve it exactly (in both directions) to generate an exact pre-classical solution. However 
the fact that approximate locally pre-classical solutions match with pre-classical solutions to 
order r y 2 and the property of local stability guarantees that a pre-classical solution defined 
from any local data and any of the approximate pre-classical solutions constructed, will 
always be 'nearby'. For instance, if we denote by e the largest of the 5j/V2k over all the 
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neighborhoods and all partitions, then we can expect any of the approximations or any of 
the exact pre-classical solutions to satisfy pairwise 

\s m -s' m \ < e V m G [M min , M max ] (30) 

Thus, our reformulated definition of pre-classicality, allows us to construct approximate 
pre-classical solutions to within e-tolerance. 
Several remarks are in order. 

(1) The requirement of admissibility of a pre-classical solution and of local stability are 
both conceptually and logically independent. It is the requirement of admissibility of a 
pre-classical solution that implies the local constancy of the coefficients (A* -2 ))™ together 
with the condition (|18|) on these constants. These are explicitly picked out by the require- 
ment of getting an approximating differential equation (Wheeler-DeWitt equation) which is 
independent of 7, m explicitly. As a by-product, this requirement provides us with locally 
approximated evolutions and also approximate, locally pre-classical solutions which are pa- 
rameterized by two free parameters. The local stability on the other hand is concerned with 
the construction of an approximate pre-classical solution by 'joining' together the various 
local approximations. Because the order of the fundamental evolution equation is always 
larger than two in loop quantum cosmology, the extension of local solutions introduces per- 
turbations of local data which only stay small if the approximated evolution is stable. By 
contrast, a similarly locally approximated differential equation does not need stability for an 
extension since the initial data are local and we have exactly the same number of parameters 
as needed for matching. 

(2) Although local stability is needed for approximate pre-classical solutions, the con- 
dition itself refers to the characteristic roots which are independent of the pre-classical 
solution. The stability property of the approximated evolution generically ensures that, 
compared to a given solution, nearby 'initial' data evolve into nearby 'final' data. If any of 
the non-pre-classical solution can also be approximated locally, then the extensions of its 
local approximations will involve similar perturbations of initial data. Note however that 
non-pre-classical solutions would generically be parameterized by more parameters (2k — 2) 
which are available for performing an exact matching and stability will not be a concern 
there (and not be necessary from a physical point of view). 

(3) The implication of stability for the roots is independent of their multiplicities (less 
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than or equal to 2k — 2) which hardly affect the growth of deviations. The dangerous 
exponential growth is prevented by stability. Strictly speaking, the definition of stability 
allows roots to have absolute values to be 1 + o((mj) _1 ). With finite domains for local 
neighborhoods, one cannot hope to derive strict equalities. However, even without strict 
equalities, stability requirement puts severe constraints on the coefficients. 

(4) It is important to note that the need for stability arises out of matching local ap- 
proximations in the adjacent neighborhoods which always introduces differences between 
extended initial data and local initial data. It only matters that the local approximations 
are obtained as solutions of an initial value problem (i.e. evolution equations). Whether the 
evolution is via a homogeneous difference equation or via a non-homogeneous one, at best 
controls the difference between the two initial data sets. For the evolution of differences 
of solutions, it is only the homogeneous part of the evolution that matters. This obser- 
vation will be particularly relevant for fundamental equations which are partial difference 
equations. 

3. Further Implications for Loop Quantum Cosmology 

Local stability provides a strict selection criterion which allows us to reduce the number 
of choices for a quantization of the Hamiltonian constraint in homogeneous models. This is 
because for a random polynomial with coefficients satisfying (114)1 it is very unlikely that all 
roots will have norm one. Also a quantization of the constraint using holonomy operators 
in some form will violate this condition in most cases. This is in particular helpful because 
models with non-zero intrinsic curvature require an additional input which is not present in 
the full theory and thus are less close to full quantum gravity. Only the Bianchi I model 
with vanishing intrinsic curvature does not require additional steps, and its quantization 
which has considerably less freedom turns out to be locally stable automatically. 

Local stability was motivated by semiclassical considerations, and it has to be fulfilled 
only in semiclassical regimes. In the Planck regime it can be violated, and it usually is in 
the presence of matter or other effects which lead to large curvature. In such a situation it 
is expected that contributions to the wave function with Planck scale oscillations become 
dominant, leading to large differences between the discrete and a continuum formulation. 
Usually, an unbounded growth of the strongly oscillating solutions when one approaches 
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the classical singularity is suppressed by a consistency condition which follows from the 
constraint equation j(| as can, e.g., be seen from the discussion in jl^ for the isotropic 
model. 

To summarize our results for the case of isotropic loop quantum cosmology: Our require- 
ments on the fundamental difference equation (|13|) are (a) the equation admits at least one 
pre-classical solution and (b) it is locally stable. The former leads to the conditions (|18|) 
while the latter requires that all the roots of the characteristic polynomial of the locally 
approximated difference equation with constant coefficients have unit absolute value. The 
conditions (j!8j) also imply that among these roots is the root A = 1 + 0(7) with multiplicity 
two and this in turn is responsible for the two parameter family of approximate, locally 
pre-classical solutions. The fact that the multiplicity is two can directly be traced back to 
the second order pole in 7 of the coefficients. A generalization to a different pole structure 
is straightforward. 



IV. PARTIAL DIFFERENCE EQUATIONS 

So far we concentrated on difference equations which were ordinary, homogeneous dif- 
ference equations. Very often, one encounters partial difference equations as in the case of 
anisotropic, homogeneous loop quantum cosmology. The rationale and the logic for contin- 
uum approximations, pre-classicality etc remains the same but the details of course differ. 
We will concentrate on the differences one encounters and keep the discussion general. The 
concrete case of homogeneous loop quantum cosmology will be used as an example; it is 
discussed in more detail in P|. 

A. Continuum Approximation 

The fundamental partial difference equation is of the form, 

00 

£<f (7)* {m ^ } = with 4;;>( 7 ) = £ iKa^ } (31) 

{ii} i=-i 
The locally slowly varying sequence is given in terms of an interpolating, locally slowly 

varying function ip(jp) by (recall j/(m/) = jirrij), 

S {mi+5inj} S{mj} 

+ 7£^5m J (9j^(p / (m / )) + -'y 2 e 2 ^25m J 5m K djd K ilj{p I {m I )) + • -(-32) 

J J,K 
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The 7-independent equations obtained from the fundamental equation become, 

E(^ H) )£} } ) v>(pV/)) = o, 



(33) 



E^ H+1) )|r;; } v^m) +^E Em h )!5 } = o (34) 

and so on. 

Once again one has to appeal to recovery of the correct continuum equation (analogous 
to the classical limit) to select the leading order partial differential equation. In the case of 
general homogeneous loop quantum cosmology (1 = 2 and no 7 _1 term), this requires the 
differential equation to be second order. 

For a non-zero solution, the bracket in eq.()33p must be zero. Using (rotational) symmetry 
conditions for the coefficients, we can reformulate this condition as 



$>(-)){-> = 



(35) 



(summing only over a single ii) for all / and all choices of ij for J ^ I. 

The second equation (jHljl either gives us a first order differential equation or leads to 
further conditions on the coefficients of the difference equation and so on. Let us assume 
for definiteness that the desired continuum equation is of second order. Then the additional 
conditions for the coefficients are, 



{ii} 







and 







(36) 



(37) 



for all J (this cannot be split into several conditions by using symmetry). 

For homogeneous loop quantum cosmology, the conditions on the coefficients are 







{*/} 



(38) 
(39) 



ij {il}l=£J} 
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and are in fact satisfied for (|1U|). In this case, the approximating partial differential equation 
is the Wheeler-DeWitt equation for the model. 

The approximating partial differential equation becomes, 

fYl ( E^(^ H) )&; } ) %k W)+*>Y, ( ) ^(/K)) 

E^ H+2) )£; } )^^)=°- ( 4 °) 

Note that the bracket in the first term can be zero for some values of J, K but it has to be 
non-zero for those J, K which will reproduce the correct continuum equation. The desired 
continuum approximation thus gives us the condition ETT j) . 



B. Local Stability 

The next step is to work with the difference equation in a local neighborhood and charac- 
terize the local continuum solutions (analogues of pre-classical solutions). These are Taylor 
expandable in 7 and lead to a hierarchy of partial difference equations. The S/^j satis- 
fies a partial difference equation with constant coefficients while the higher order terms in 
the Taylor expansion are determined recursively from non-homogeneous partial difference 
equations with constant coefficients. In a local neighborhood, we can take a short cut and 
directly take an ansatz for a local pre-classical solution as, 

(0) (K))/) +i^5m I b ( i\(mv) I ) + 1 2 ^5m I 5m J c i ?}((m ) I ) + --- . (41) 
/ i,j 

Substitution then implies that ((m ) 1) , b^p ((m ) 1) are free parameters analogous to 
the case of ordinary difference equation. Now however there are additional possible free 
parameters, for example the cfj((m )j) for those J, J for which the first bracket in eqn ()40|) 
is zero. Similar argument implies that there could be free parameters at higher orders in the 
Taylor expansion (|41j) . This is different from the case of ordinary difference equations where 
the higher order coefficients are all determined in terms of the a^°\ b^. The local continuum 
solution is thus parameterized by potentially infinitely many free parameters. Extension of 
a local continuum solution to the adjacent 'cells' of the given local neighborhood will be 
more complicated. 
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In the previous section, for formulation and analysis of 'local stability' we concentrated 
on the locally approximated evolution, equation ([27)1 . The analogue of this equation is the 
equation satisfied by namely, 

£(^- 2) )^4t + M = - ( 42 ) 

{</} 

The (mo)/ dependence of the coefficients is suppressed. 

Unlike ordinary difference equations, this is not an evolution equation and we do not 
have a handle on its general solution. We can think of fixing mj + ij,I ^ 1 and allowing 
variation only along the 'first direction' and see if an evolution equation can be arrived at. 

Let the lattice be D- dimensional. Let us introduce the notation: 

Z%' 1 l3 '-' tD (m 2 ,m 3 ,---,m D ) := sS, m2+i2 ,..., mj3 -B , 

^ 'M S {-mi+i I } 

h 

= £(^ ( - 2) W,te Z^r^(m 2 ,m 3r ..,m D ) (43) 

h 

The equation ()42|) then can be expressed as, 

X i 2 ,is,-,in = o. (44) 



E 

12:13, — ,in 



We can separate a particular set of values 1%, i^, • • • , i* D and view the above equation 
as an equation for X l2 ' l $' '"' 1 d. Provided, X l2,l3, "' ,lD , involves a non-trivial combination of 
^mf+ii ' lD ( m 2, ■ • • > m D) (he. more than one mi + ii occurring in the sum), we can think 
of the equation as an ordinary, non- homogeneous difference equation with respect to m 1; 
for the -Z^i+'u %D { m 2-, "13, • • • , mc). Thus we can interpret the partial difference equation 
as a (non-homogeneous) evolution equation for evolution along a line, parallel to the m\ 
direction and passing through m 2 + i\, m 3 + i\ y • ■ • , mr, + i* D in the lattice. 

For example, in the context of homogeneous loop quantum cosmology, the lattice is three 
dimensional and the indices ij range from —2 to +2. The locally approximated partial 
difference equation (}4*2*|) turns out to be of the form, 



7/ W JO) I ( ) _ ( ) _ ( ) 

U[mi) I S mi m2+ 2,m 3 +2 ' ■ S mi,m 2 -2,m 3 -2 S mi,m 2 +2,m 3 -2 S m 1 ,m 2 -2,m 3 +2 

_L_ At \( (°) (°) _ (°) _ (°) 

' a \ m 2) I S mi+2,m 2 ,m 3 +2 ' S m 1 -2,m 2 ,m 3 -2 S m 1 -2,m 2 ,m 3 +2 S m 1 +2,m 2 ,m 3 -2 

+ rff m ^ + s (°) _ s (°) _ s (°) ) - (A5) 

T uyiibz) I S mi+ 2,m 2 +2 ) m3 ' a mi-2,m 2 -2,m 3 s mi-2,m 2 +2,m3 s mi+2,m 2 -2,m3 ) ~ ' * ' 
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where d(m) m m _1 for large m. 

The equation is invariant under permutation of the three directions, so it is enough to 
consider 'evolution' along the 'mi '-direction. For the (^2,^3) set we have the eight pos- 
sibilities (2, 2), (2, -2), (-2, 2), (-2, -2), (2, 0), (-2, 0), (0, 2), (0, -2) and the corresponding 
^m'+ii- The corresponding X' 2 ' 13 involve non-trivial combinations of Z J s only for the four 
sets (2,0), (-2,0), (0,2), (0,-2). For example, X 2 ' = d(m 3 ) (Z^ +2 -Z^°_ t ) is a non-trivial 
combination while X 2 ' 2 = (f(mi)Z^ is a trivial combination. Clearly, the partial difference 
equation can be viewed as an ordinary, non-homogeneous difference equation for any of the 
four non-trivial combinations. For each of these, the order of the evolution equation with 
respect to m\ is 4. This illustrates how the locally approximated partial difference can be 
viewed as (different) non- homogeneous evolution equations for different combinations, X's. 

Matching various local solutions across adjacent cells requires, in particular, matching 
along various directions connecting the adjacent cells. Now the arguments used for the case 
of ordinary difference equations can be invoked and issue of local stability becomes relevant. 
Note that these arguments are relevant only for the non-trivial combinations. For trivial 
combinations, local solutions in adjacent cells, across a cell boundary are decoupled and 
hence matching requires only equating the values on the cell boundary. 

As observed before, the criteria for local stability depend only on the equation being an 
evolution equation and not on its (non-)homogeneity property. The criterion in terms of the 
characteristic roots of the associated homogeneous equation, X' 2 '^'' " ,%D = 0, then remains 
the same. This criterion is of course to be applied to all the evolution equations for different 
sets of («2, 13, ■ ■ ■ 1 Id) and along each of the directions. These details depend very much on 

(—2) 

the precise coefficients i2 ... • and are to be analyzed on a case by case basis. 

For the equation of homogeneous loop quantum cosmology, eq. (j43j) . the characteristic 
roots of the evolution equations are ±1, ±z which all have absolute value 1 thus satisfying 
the criterion of local stability jsj . 

To summarize: The rationale for continuum approximation and local stability remains 
the same as in the case of ordinary difference equation. In particular, with regards to local 
stability, one still deals with the locally approximated partial difference equation. Whenever 
this equation can be viewed as an evolution equation, we require this evolution to be locally 
stable as defined in the subsection IHI Bl Since there could be multiple ways of viewing the 
partial difference equation as an evolution equation, local stability of all these evolutions is 
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required. 



V. CONCLUSIONS 

In this paper we have taken the view point that the fundamental equations of a the- 
ory arise from a discrete formulation and are thus difference equations. Typically such a 
fundamental discrete formulation arises from a quantization procedure which is not free of 
inherent ambiguities. A basic requirement on any such discrete formulation is that a contin- 
uum approximation should be available in suitable regimes. The parameter(s) controlling 
the discrete structure, being physical, cannot be taken to zero in arriving at a continuum 
approximation. A continuum approximation is then to be derived by invoking a suitable 
notion of pre-classicality without using unphysical, mathematical limits. In the context 
of isotropic loop quantum cosmology, we have shown how the definitions of slowly vary- 
ing solutions together with 7-independence of the corresponding approximating differential 
equation suffices. Furthermore the formulation is powerful enough to put constraints on the 
coefficients of the fundamental equations as in equation ()18j) . These are however necessary 
conditions for a continuum approximation to be admissible. The sufficiency can also be 
demonstrated in a well defined approximate sense by a constructive procedure provided, in 
addition, the condition of local stability (J29|) is satisfied. 

More generally, the fundamental difference equation will be a partial difference equation. 
These arise for instance, in the context of the more general homogeneous models. The logic 
of admissibility of a continuum approximation is still valid and is seen to lead to similar 
conditions on the coefficients. One is then lead to a locally approximated partial difference 
equation. For the favorable case when the locally approximated partial difference equation 
can be viewed as a non-homogeneous, evolution equations along some direction for some 
1, one can extend local solutions across neighboring cells in a manner similar to the 
case of ordinary difference equation. The issue of stability of local approximations arises 
again for exactly the same reasons and leads to the definition of local stability exactly as 
before. The local stability criteria however are to be applied to all the evolution equations 
implied by the locally approximated partial difference equation. For the fundamental equa- 
tion of the diagonal, homogeneous models, these evolution equations turn out to be stable 

G 
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When the fundamental difference equation is derived from a (quantum) theory, there are 
usually free parameters or other choices (e.g., quantization ambiguities). In such a case 
the conditions derived here can serve as powerful restriction criteria which go beyond the 
usual requirement that the quantization have the correct classical limit. Thus, they provide 
consistency checks for possible quantization approaches and restrict possible choices within 
a given approach. In this light, it is quite remarkable that the loop quantization of the 
homogeneous Bianchi I model, which is very close to the quantization of the full theory and 
does not allow many choices, turns out to satisfy our conditions automatically. This can be 
regarded as a positive consistency check for the full theory. For other homogeneous models, 
which have non-zero intrinsic curvature and are therefore not as close to the full theory, 
there are more choices which in particular allow violations of the conditions. In this case the 
conditions serve as strong selection criteria since it is possible to find quantizations fulfilling 
them for all those models [31. 
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APPENDIX: AN EXAMPLE OF AN UNSTABLE QUANTIZATION 

In this appendix we will focus on the isotropic subclass of the Bianchi-I and Bianchi-IX 
models and examine the stability properties of the evolution equations resulting from the 
quantization of the constraint operator as given originally in 9]. Again for our purposes, it 
is sufficient to consider only the vacuum case without the cosmological constant. 

Let us recall some basic expressions: An orthonormal basis of states in the connection 
representation is given by (c\n) analogously to (J7J), which also forms an eigenbasis of the vol- 
ume operator. In the isotropic case the convention V\ n \-i = (§7^p)^ \/(M ~~ 1) \ n \ (\ n \ + 1) 
has been used. As in the full theory jlfij . the constraint H = —He + P consists of two 
parts which can be quantized separately (in the main text we quantized the full constraint 
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directly which cannot be done in the full theory): 

H E = AiijKil)- 1 e'^trihjhjhfhfh^hKlhx 1 ^}) 

I,J,K 

where h[i t j] := Y[k(^ k ) Cij with the structure constants Cf 3 = 6uk for the closed isotropic 
model (Bianchi IX) and zero for the flat model (Bianchi I). The second part of the constraint 
operator is obtained from the first one by 

P = - 8 z(i + 7 - 2 ) /t - 1 ( 7 4)- 3 t IJK ti{h I [h-\k]hj[hf,k}h K [h- K \v}) 

I,J,K 

with K = -| 7 - 2 ^- 1 [ J ff E) V] sgn^deTa} ) . 

Upon transforming to the triad representation by \s) := ^2 ne z^n\n), the constraint 
equation H\s) =0 gets translated to a difference equation for the coefficients S n which are 
just complex numbers in the vacuum case. The order of this equation is 16 for Bianchi-I and 
20 for Bianchi-IX. Since the coefficient of So is always zero in the factor ordering chosen for 
the constraint operator, one can ignore the sgn factors which enter through the definition of 
the K operator i.e. equivalently absorb them in the S n 's. Explicitly the equations are: 

For Bianchi-I: 




For Bianchi-IX: 




1 -)- 

+ ( J ^2 ( ykl=£M — V\n-2k\-2 j A n ^ k 2k S n ^2k (A-2) 

fc=-5 

The definitions of and A l n are given in 

Observe that we must have infinitely many S n to be non-zero for a non-trivial solution 
and the number of non-trivial independent solutions is reduced by one in both cases since 
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the coefficient of Sq is zero (see p|). We are interested in the admissibility of a continuum 
approximation and in particular in checking if the evolution is locally stable. 

To explore these properties, we consider the large n behavior of the coefficients, S n . For 
defmiteness we take n to be positive and large which gets rid of the absolute values in the 
eigenvalues of the volume operator. Introduce the notation: 



a\n) := ^ (A.3) 
o 



For large n one gets, 



V*p - a ( 1 + 2^T + 8n» + '") (A ' 4) 

ftf — l + o(±) (A.5) 

For Bianchi-IX one needs A l n defined in terms of := [K. q n+1 — JC q n _ 1 )/2 for q = 
±1,±3,±5. These become, 

^ - =F (^) { =F »» + I =F » + .(!)} (A.6) 

- ±(-3g-){ t y» + y t y» + »(-)} < A -7) 
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l±i\ r 45 225 45, /1\-| 

=F — 7i H =F — b + o [ — )> (A. 

2 8 2 Vn/J V 



— + °(r) (A.9) 



4 / Vra 

*a - -(^) + (a.io) 



5(1 ±ih /l 



With these, the difference equations become: 

= S n+4: - 2S n + Sn-4 (S n+8 - 2S n + S n - S ) (Bianchi-I) (A. 12) 

= (1 - i)S n+5 + (1 + i)S n - 5 + (1 + 5i)S n+3 + (1 - 5i)S n - 3 - 2(1 + i)S n+ i 

1 + 7 ~ 

fc=-5 

The coefficients A™ satisfy the unitarity conditions A~ m = (A™)* and are obtained as: 



i,-2 s 

-2(1 - i)S n „t + (—^-) ( Yl A n-2k S n-2k) (Bianchi-IX) (A.13) 



A™ = +—i 
32 

8 45 15 

K ~ "16 ~T 
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A n = -T + T» 

n 4 4 

< = +f (A.14) 

Note that in the large-n limit, the equation has reduced to a difference equation with 
constant coefficients and our first necessary condition for admissibility of continuum approx- 
imation is satisfied. The condition of Eq. (fT5|) however is satisfied by the Bianchi-I case but 
not by the Bianchi-IX case. Therefore we will not get the Wheeler-DeWitt equation in the 
leading approximation for the latter. 

A simple check on local stability is to obtain the characteristic roots to the leading order 
in 7. This picks out the polynomial obtained from the coefficients in the term multiplying 
T 2 . 

In both cases of Bianchi types, the corresponding characteristic polynomial is a polyno- 
mial in z 2 . Hence both ±z are roots. The coefficients of the polynomial also come in complex 
conjugate pairs such that if z is a root then so is 1/z*. Furthermore, these polynomials are 
perfect squares, so each root is a double root. Corresponding to each such root, there will 
be solutions which will behave as z n (an + b). 

For Bianchi-I, in the continuum approximation, the roots of the characteristic polynomial 
are the 8 th roots of unity each being a double root. Clearly all those roots have absolute 
values equal to 1 and exactly one root equals 1. Our condition of local stability is thus 
satisfied (see also for the continuum limit). 

For Bianchi-IX on the other hand, the roots are explicitly obtained as: 



A = 0.24077449441476 - i 0.47063359530671 
Ai = 0.29768070243356 + i 0.95466549083889 
A 2 = 0.8 -i 0.6 

A 3 = 0.86154480315169 -i 1.68403189553218 
A 4 = -1.0 



|A | = 0.52864765032298 
[Ai| = 1.0 
|A 2 | = 1.0 

|A 3 | = 1.89161911414729 
|A 4 | = 1.0 (A.15) 

The actual 10 double roots are zf = ±a/A7. Observe that 1 is not a root of the charac- 
teristic polynomial. Hence the approximating differential equation equation will also have 
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lower derivative terms. Moreover we do have roots whose absolute values are not equal to 
one, thus violating local stability. 

Thus, the original quantization in jjj for the closed isotropic model is not admissible by 
our criteria. In 12], a different quantization of the closed model has been derived with the 

n 

methods of |3J, providing an admissible quantization for this model. 
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